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Abstract 

> . 

■^j- , According to one-loop perturbation theory, fermions whose masses 

are totally generated from Yukawa couplings do not decouple in the 
^ ■ heavy mass limit. We investigate this issue nonperturbatively in the 

strong coupling regime of the chiral U(l) mirror-fermion model in four 
dimensions. Our numerical results, obtained on 6 3 • 16, 6 3 • 24 and 8 3 • 16 
lattices, indicate nondecoupling of heavy fermion and mirror-fermion, 
thus supporting the one-loop picture. 



<D . 

^ 1 Introduction 

> : 

The Higgs and heavy fermion sectors in the Standard Model have received a lot of at- 
tention from the lattice community recently. The main difficulty in numerical studies 
of them is how to deal with chiral fermions on the lattice. Different approaches were 
proposed M. The approach adopted in this letter is the mirror-fermion method of 
Montvay @. 

In the symmetric (PM) phase of the mirror-fermion models, the fermion and its 
mirror partner are always degenerate. Only in the broken (FM) phase where the 
scalar field develops a nonzero vacuum expectation value (VEV) can one manage to 
decouple the mirror-fermion. 

One way to decouple the mirror-fermion is to decouple it like a right-handed 
neutrino |J. This is possible due to the shift symmetry of the action at vanishing 
mirror-fermion Yukawa coupling. Another possibility is to give the mirror-fermion a 
very heavy mass (around the cut-off scale or almost as heavy as the lowest doubler) 
and hopefully they will have vanishing effects on the scalar field in this limit ||]. 
This kind of decoupling is based on the decoupling theorem || saying that when 
the particle has a mass much higher than the physical scale, it will have very small 
influence (e.g. radiative corrections) on the "physical world". However, it can be 
easily seen at one-loop level that the contribution of the particle to renormalized 
quantities of other particles will not be suppressed by its huge mass if its mass is 
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generated from the Yukawa coupling, i.e. if its mass is generated through the mech- 
anism of spontaneous symmetry breaking (SSB). This is the so-called phenomenon 
of nondecoupling. 

Nondecoupling of heavy fermions in theories with SSB has been discussed in 
several papers H [7], |8|]. So far, all arguments are within (one- loop) perturbation 
theory (except in |J where large-N expansion is used). Nondecoupling beyond one- 
loop is still not clear. We decide to study this issue in a nonperturbative way in our 
£7(1) scalar-fermion model with explicit mirror-fermions. We still ignore the gauge 
field assuming that this approximation will not change the picture qualitatively. 
Hence, the scalar field is our "physical world" . 

2 Action and Renormalized Quantities 

The Euclidean lattice action for the mirror- fermion model with chiral U(l) symmetry 
in four dimensions in the usual normalization convention for numerical simulations 



is 



r 1 2 r 

s = Ysi^t^x) + A [(<pi<p x ) - l\ - /cj](^«) + ^ x [{xM + {^ x Xx) 

X fjj 

-KY^lfflx+pjM + {Xx+nlnXx) + r ({Xx+n^x) ~ {Xx^x) + (JP x+fl Xx) ~ @ X X. 
p- 

+ G f [lp x ((j>lx - 17502*)^*] + G X lXx(<Plx + ilsMXx}} • (1) 

Here K is the fermion hopping parameter, r the Wilson-parameter, which will be 
fixed to r = 1 in the numerical simulations. The complex scalar field is <p x = <pi x + 
i(f>2x, and ^ x = (ip x , Xx) stands for the mirror pair of fermion fields (usually ip is the 
fermion and x the mirror-fermion) , yj M sums over eight directions and 7_ M = — 7 M . 
In this normalization the fermion-mirror-fermion mixing mass is \i^, x = 1 — 8rK. 

Renormalized quantities in the FM phase are defined as follows: 
For a given configuration of the scalar field, we measure 

1 5>. = Me fa . (2) 



L 3 T 



X 



The bare VEV is defined as 

v = M). (3) 

For each given configuration, the longitudinal (<f>Lx) and transverse (4>tx) scalar field 
components are defined as: 

<f>' x = e~ ia (p x = (pLx + i <Ptx ■ (4) 
At the same time, the fermion fields should be rotated: 

V> ' Lx = Z~ la ^L X , V'ito = ^Rx , $'Lx = ^L X e ta , i>' Rx = IpRx , 

Xrx = e~ ia XRx , x'lx = Xl x , x'r x = XRxe ia , x' Lx = Xlx ■ (5) 
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The reason why we make transformations on all fields for each configuration 
is because there is actually no SSB on a finite lattice. Numerical measurements 
of non-symmetric quantities on a finite lattice will always end up with zero. We 
therefore need to make the above rotations and define the bare VEV as the length 
of (ft in eq.(2). This quantity will converge to the bare VEV of the scalar field in the 
thermodynamic limit ||. In the same limit, <pL X — v and <pTx wm converge to the 
Higgs and Goldstone fields. 

In our numerical simulation in the FM phase, renormalized quantities for the 
scalar field are defined by the timeslices of the longitudinal and transverse correla- 
tions ^ 

S Lt = 71 XX(<?W0LOo})c = 71 ^2{(4>Lxt<t>LO0) ~ V 2 ) , 
^ x L x 

S Tt = 71 XX(0Txi</>TOo))c = 71 XX(<^xt</>TOo)) • (6) 

Here it is taken into account that (4>l,x) — v an d (4>Tx) — 0. 
The physical Higgs boson mass is obtained from 

= a + 6(e- mi * + e - mL{T - t] ) . (7) 

The constant a is due to finite size effects [flQ[1 . 

The wave-function renormalization constants and Zt are defined as 

m| + 4sin 2 ^) ^ cos(^-t) (<Plx^lo)c , 
Z T =A sin 2 ^ X! ^(S^) (^0to) • (8) 

X 

Because of the existence of the massless Goldstone boson in the FM phase, longi- 
tudinal quantities, e.g. Zl, cannot be defined at zero momentum. Although the 
transverse mode does not suffer from infrared problems, Zt cannot be defined at 
zero momentum either because <pTx corresponds to a massless particle which has no 
rest frame. 

The renormalized VEV, denoted by vr, is defined as 

VR = -7=- • (9) 



Renormalized fermionic quantities can be defined from the 2-point function (p) 
of the transformed fermion field ty' = (■?//, x')' 

A*(p) = e-*"^ (%%) . (10) 

X 

At momenta p = (0, 0, 0,p^) where p4 is small, we have 

A*(p) ~ A-i sin(p 4 ) -IaB , f*(p) = A^\p) ~ M + i sin(p 4 ) "74^ , (H) 
with 

M = [A + sin 2 ^)^^" 1 ^) -1 , N = A~ X BM = MBA' 1 . (12) 
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1/2 

The fermion wave-function renormalization matrix has to transform N to 
the unit matrix: 

(4 /2 fiv4 /2 = i (13) 

where superscript T means transpose of the matrix. The renormalized fermion mass 
matrix is obtained as 

M R ee (4 /2 ) T M4 /2 = ( ^ ) , - ^ • (14) 

\ Mii Vrx J vr 

In general, //r / 0, we therefore need to diagonalize the renormalized mass 
matrix Mr to obtain the physical fermion states , which will be mixtures of the 
original ip and x states with the mixing angle ocr defined later. If the diagonalized 
mass matrix is denoted by M p h ys , we have 



■Mphyg 



fJ>W 
H2R 



Hir = Hity cos 2 (a R ) + ii R sin(2a R ) + \i Rx sm 2 (a R ) ee G 1R v r , 
H2R = Hr x cos 2 (a R ) - ii R sxn(2a R ) + /Ir^ sm 2 (a R ) = G 2 rv r , 

2 afl = tan- 1 ( ^ R ) . (15) 



Our convention is always that after diagonalization of Mr, the second physical 
fermion state is the heavier one, i.e. \(j<2r\ > l^ifll- The sign of the fermion mass is 



insignificant, because it can be changed by an appropriate 75-transformation [p"T|j . 

The 3-point renormalized Yukawa couplings coupled to the transverse scalar field 
(pT (i.e. the Goldstone boson)[] are defined by 



1 ^G ( rI \ A _ jLl\.f,..\ 7 ->r2 r , tl , l7 ->r2,T 

-ilsG%l 



6 k ,„ p+q = ^fR(p 4 )Z- 1/2 GV(Z- 1/2 ) T f R (q 4 ) , (16) 



where k±, p±, are the 4th components of the momenta of Goldstone boson, fermion 
and anti- fermion, respectively, and fc| is 4sin 2 (k 4 /2). We have set the spatial com- 
ponents of all momenta to zero. The appearance of the Kronecker-delta above is 
due to energy-momentum conservation where 

G (c) = _j_ J2 e- lkiXi e~ ipAVA e iqAZA (^ Tx ^' y K 

x,y,z 



is the connected part of the 4>t-^'-^' 3-point Green's function. Since {4>tx) = 0, 
the above connected Green's function is equal to the disconnected one. In our 
simulations on L 3 • T lattices we choose 

2ll 7T 7T 

k i = — , Pi = -- , q± = - . 



1 Here one can think of the Goldstone field as the longitudinal component of the gauge field 
once the system is gauged. 
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After carrying out all the matrix multiplications, we get G^ and G#*. The expres- 
sions are too voluminous to be displayed here. 

If fiR ^ 0, one more diagonalization is necessary to obtain the 3-point renor- 
malized Yukawa couplings for the physical fermion states, denoted by G\ R and G2R 
respectively We have 

Gfl = G%cos 2 (a R ) - G^sm 2 (a R ) , G^ = G^cos 2 (a R ) - G^sin 2 (a R ) . (17) 

We then define the following ratios: 

r (3) ^(3) ^(3) „(3) 

p __ <^RiP p _ <fJ_R . p __ ^R X p _ ( f2R n a\ 

^Rij> ^1R tjRx 

Note that renormalized fermionic quantities defined from eq. (|TTD to eq. fll8|) are valid 
at the zero corner of the Brillouin zone, i.e. they represent properties of the two 
physical fermions. One can actually define corresponding quantities for fermions at 
other corners in the Brillouin zone. For instance, to obtain quantities for the fermion 
doublers at p — (eiir, e 2 7T, e 3 7r, e^ir) ; ei i2 ,3,4 = 0,1, one simply makes the following 
transformation on the rotated fermion field 

y _ ^_^yix 1 +e 2 x 2 +e 3 X3+e 4 X4 ^jl (19) 

Then fermionic renormalized quantities for the two doublers sitting at p will be 
obtained from eq.flTTD to eq.([T8|) by substituting the transformed fermion field for 
the original one. Renormalized fermionic quantities for the doublers will have the 
same notations as their corresponding quantities for the two physical fermions except 
that we attach a bracket in which we denote where in the Brillouin zone they are 
sitting. For example, ratios Ri and R2 for the doublers at one-7r corner are denoted 
as i?i(l7r) and R 2 (It!') respectively. 

As we mentioned earlier, one possible way of dealing with the mirror-fermion is 
to give it a very heavy mass and hopefully it will have vanishing effects on the scalar 
field in this limit Q. This can be easily seen at one-loop level where the scalar 
propagator receives radiative correction from the fermion loop whose contribution 
is proportional to the square of the ratio of the two Yukawa couplings (e.g. R^ 
in (|r5D). In the heavy fermion mass limit, this ratio will be suppressed by the 
factor of fermion mass. However, in a theory with SSB, the Yukawa coupling itself 
is proportional to the fermion mass 0. In this case, the ratio say, R^ will be a 
constant and the radiative correction will not be suppressed by the heavy fermion 
mass. Instead, it will be proportional to l/v R , hence the heavy fermion will not 
decouple from the real world. It is necessary to go beyond one-loop and investigate 
this issue in a completely nonperturbative way, i.e. on the lattice. 

Since fi R 7^ in general, the relevant ratio to measure for the study of decoupling 
of heavy mirror-fermion is R2, which is very close to R x when n R ~ 0. According to 
one-loop calculation, 3-point Yukawa couplings defined in ( PS|) coincide with G R s 
defined by the mass-to- VEV ratio in (|14D. When the Yukawa couplings are strong, 
they may differ from each other and G R s in ([14]) cease to be a good definition for the 
renormalized Yukawa couplings. There, the better definition should be the ones in 
(|16|). What is interesting to see is whether the two Yukawa couplings still coincide 
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with each other beyond one-loop level. If this is the case in the strong-coupling 
regime such that R2 = 1.0 within the error, then we take it as the nonperturbative 
signal for nondecoupling. 

At tree level, it is clear that lir for the doublers will have contributions from the 
Wilson term. Therefore we expect that masses for the "physical" (i.e. diagonalized) 
doublers will not be totally generated from Yukawa interactions, and decoupling 
theorem should apply to all of them. Corresponding ratios R\ and R2 are also 
measured to see whether doublers are decoupled as expected. 



3 Numerical Simulation 

The Monte Carlo simulations were performed by the unbiased Hybrid Monte Carlo 



method |12| . Therefore, the fermions have to be doubled by taking the adjoint of the 
fermion matrix for the second species. (The fermionic part in ([!]) is given for a single 
mirror pair of fermions.) The number of leapfrog steps per molecular trajectory 
was chosen randomly between 3 and 10. The step size was tuned to maintain an 
acceptance rate around 75%. The necessary inversions of the fermion matrix were 
done by the conjugate gradient iteration, until the residuum was smaller than 10~ 12 
times the length square of the input vector on the 6 3 • 16 and 6 3 • 24 lattices. On the 



8 3 • 16 lattice, a precision of 10~ 13 was necessary. 

In the simulations, we set A = 00 and G x = —0.1, —0.3, —0.6 and —1.0, while 
Gj/j is kept at 0.1. The mass parameters k and K are tuned such that the mixing 
mass in the FM phase is zero or small, and the renormalized fermion mass li\r is 
around 120 GeV in physical units, mimicking the top quark. (The scale is set by 
vr = 246 GeV.) We let and G x have opposite signs because this combination is 
preferred by the condition of reflection positivity PJ. The mixing mass lir is tuned 
to be small in order that masses of the physical fermions are basically completely 
generated from SSB, which is the case we interested in. The most difficult task in 

this numerical simulation is the tuning of K to have lir — 0.0. (Our criterion for 

(3) (3) (3) 

/jLr ~ 0.0 is that lir^, llr x , G r L and G R ' do not differ from liir, Li2R, G\r and 

G2R respectively by more than 1%. In this case, the second physical fermion state 
is basically the mirror- fermion.) When G^ ■ G x = 0, K c (the value of K at which 
lir = 0.0) is at 1/8. When the bare Yukawa couplings are small and G^ ■ G x < 0, 
K c will be slightly larger than 1/8 according to one-loop bare perturbation theory 
[Tl[| . However, we find that even when — G x > 0.6, K c still keeps growing as is 



qualitatively predicted by one- loop result. 

Most of the runs were done on the 6 3 • 16 lattice. Since it is not a very large 
lattice, we need to check the finite size effect. This can be estimated by one-loop 
renormalized perturbation theory on the lattice. Formulae of one-loop bare per- 
turbation theory on the lattice were already presented in 0]. Starting from those 
equations, one can obtain one-loop renormalized perturbation theory on the lattice 
where counterterms are introduced on the infinite lattice to eliminate ultraviolet 
divergences. For example, the bare fermion mass is given by eq.(74) in ||. By 
assuming that = 1/2K and stays constant as we change the lattice size (which is 
a good approximation to our numerical data), we obtain the finite size shift of the 
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renormalized fermion mass A/i^ = fii R (oo, oo) — fti R (L,T) at fi R = as 

1 1 



Afi 1R = (jl 1r G 2 m A / [ 2 - -] .0(g)- 1 (<f + /x 2 R ) (20) 

J q L q -+- m L q 



where 



Notice that all quantities appearing on the right-hand side of eq . (|20"D are those in 
the infinite volume limit. Similarly, one can derive one- loop results of finite volume 
shifts for other renormalized quantities. For instance, one should start from eq.(68) 
in H to derive Am^. Those formulae are too voluminous to be displayed here. 

According to our previous experience, scalar quantities usually suffer from finite 
size effects the most due to small fermion mass fi\ R . From one-loop results, we do 
find that at G x = —1.0 the effect amounts up to 10% for mi, while G Rx and G Rx 
are smaller by about 5% in the infinite volume limit. 

Also, it is important to have some idea about the finite cutoff effect since some- 
times the cutoff is low. The way we estimate finite cutoff effect is to compare 
one-loop /3-functions in the continuum and on the lattice. One-loop /3-functions in 



the continuum are given by eqs.(83), (84) in |pT| . Starting from eqs.(74) to (79) 
in one can obtain one-loop formulae for Gr's on the lattice. By differentiation 
with respect to the energy scale, we will obtain their one-loop lattice /3-functions. 
The derivation and final results are very lengthy and will not be stated here. Those 
one-loop results show that at G x = —1.0, as we go to smaller mass scale say, from 

(3) 

nriL = 1.3 to 1.0, both G Rx and G R ' decrease. But the decrease of them is faster for 
the continuum /^-function. We estimate that at mi = 1.0, the cutoff dependence of 
both Gr x and G Rx is around 10%. Of course, this is just a crude estimate. We also 
see that according to one-loop calculations, finite volume and cutoff effects for G Rx 
are very close to those for G R ' (with the same sign). This indicates that the ratio 
i?2 should be quite stable against those effects at one-loop level. 

It is always better to check the finite volume and cutoff effects numerically at the 
same time, since one-loop formulae break down in the strong-coupling regime. Due 
to limitations on CPU time, the numerical check was carried out only at G x = —1.0, 
which is the most relevant point, on 6 3 • 24 and 8 3 • 16 lattices. There we had to 
go to the negative-K region in order to increase the cutoff. Although reflection 
positivity cannot be proven when k < ||, data obtained there did not show 
anything abnormal. We simply assume that reflection positivity is still preserved at 
least in the region with small negative k. 



4 Conclusions and Discussions 

Our numerical data are presented in tables 1, 2 and 3. From data at points d, D, e, 
and E, finite size effect for at k = —0.01, G x = —1.0 is quite strong as expected. 
It is even stronger than predicted by one-loop formula. The relevant quantities G Rx 
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(3) 

and G R ' show smaller finite volume shifts, and R2 is quite stable as we increase the 

(3) 

lattice size. Finite cutoff effect in G R ' is smaller than predicted by one- loop result. 
But fluctuations at point E are quite large. 

Readers should note that at point cc, the mixing mass hr 7^ 0, which means 
the mass of the second physical fermion has contribution from the ordinary (off- 
diagonal) mass term. Therefore, it is not surprising that R2 < 1.0 there, showing 
"partial" decoupling. For all other points, parameter K is tuned such that /jlr is 
quite close to zero. There the ratio R2 is always equal to 1.0 within the error up to 
G x = —1.0 which, in our model, is a strong coupling case. Hence, it is obvious from 
our data that the heavy mirror-fermion does have its renormalized Yukawa coupling 
proportional to the mass and will not decouple from the transverse scalar field. The 
one-loop picture survives the strong Yukawa coupling limit. The idea of decoupling 
the mirror partners by giving them large masses does not work. Since the action 
in (|l]) has # G x , 75 # —75 symmetry, it is obvious that the heavy fermion 
itself will not decouple either. It is also obvious that doublers at one-7r corner do 
decouple. Fluctuations of data of higher doublers are too large. But we believe that 
all doublers decouple like the lowest ones do. 

We think that these conclusions also apply to the SU(2) version of our mirror- 
fermion model because its qualitative behaviour appears to be similar to that of the 
£7(1) model |TB| . We conjecture that in other scalar-fermion models, heavy fermions 



do not decouple either. 

A direct measurement of S defined defined in eq.(2) in || will be overwhelmed 
by fluctuations for our models. Reasonable statistics can be reached if we study 
a scalar-fermion model with fewer internal degrees of freedom like the model with 
staggered fermion coupled to one-component scalar field. 

Once the chiral gauge field is put in, the shift symmetry at G x = fi = in 
our models is broken. This means that the mirror-fermion can no longer be de- 
coupled as right-handed neutrino in the presence of gauge field. Some other means 
for decoupling should be used. We conjecture that the mirror-fermion will remain 
coupled to the scalar field in the heavy mass limit when the system is gauged. The 
more serious problem there is actually the renormalized gauge coupling. In order 
to preserve the chiral gauge symmetry in the action, both ip^ and xr have to be 
gauged. It is possible that in the FM phase once ipL is coupled to the gauge field 
(i.e.: the real world), \R is & lso coupled at the same time as is shown at one- loop 
level. If this is true nonperturbatively, then mirror-fermion cannot be decoupled 
eventually and the gauged mirror-fermion models have trouble in reproducing the 
minimal Standard Model in the continuum limit. This issue will be explored in 
details in the future. 

I would like to thank I. Montvay and G. Minister for useful discussions. 
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Table 1: Data on renormalized VEV and masses are presented. Points with lower- 
case letters are obtained on 6 3 • 16 lattice while point D is on 6 3 • 24 lattice, point 
E on 8 3 • 16 lattice. Each point has about 4000 trajectories for equilibration, and 
around 10000 to 12000 for measurements. 









K 


vr 


m L 


Hr 


UlR 


\^2R\ 


a 


-0.1 


0.160 


0.1256 


0.330(11) 


1.06(4) 


0.033(6) 


0.2008(4) 


0.2013(4) 


b 


-0.3 


0.126 


0.128 


0.205(11) 


0.77(6) 


0.041(10) 


0.118(2) 


0.343(6) 


c 


-0.6 


0.083 


0.131 


0.28(2) 


1.01(6) 


0.040(12) 


0.1426(9) 


0.804(8) 


cc 


-0.6 


0.083 


0.1256 


0.214(17) 


0.74(6) 


0.226(10) 


0.183(1) 


0.729(18) 


d 


-1.0 


0.0 


0.1345 


0.248(3) 


1.31(8) 


0.052(20) 


0.143(1) 


1.29(3) 


D 


-1.0 


0.0 


0.1345 


0.240(8) 


1.21(7) 


0.037(37) 


0.144(8) 


1.22(7) 


e 


-1.0 


-0.01 


0.1345 


0.229(31) 


1.35(10) 


0.070(31) 


0.124(2) 


1.18(3) 


E 


-1.0 


-0.01 


0.1345 


0.231(20) 


1.03(17) 


0.066(20) 


0.118(2) 


1.026(42) 



Table 2: Data on renormalized Yukawa couplings for the fermion and mirror- 
fermion. Labels are the same as in table 1. 





Gr4> 


r (3) 


Gr x 






R2 


a 


0.60(2) 


0.60(3) 


-0.60(2) 


-0.59(3) 


-0.59(3) 


0.99(4) 


b 


0.56(3) 


0.56(4) 


-1.66(9) 


-1.70(9) 


-1.70(9) 


1.02(3) 


c 


0.50(4) 


0.48(4) 


-2.87(21) 


-2.89(22) 


-2.89(22) 


1.01(3) 


cc 


0.57(3) 


0.59(5) 


-3.13(19) 


-3.07(19) 


-2.90(18) 


0.85(4) 


d 


0.57(5) 


0.63(6) 


-5.19(49) 


-5.09(48) 


-5.08(48) 


0.98(6) 


D 


0.60(9) 


0.57(11) 


-5.10(27) 


-5.03(32) 


-5.03(31) 


0.99(6) 


e 


0.52(4) 


0.53(5) 


-4.92(33) 


-4.67(34) 


-4.66(34) 


0.95(6) 


E 


0.51(4) 


0.525(54) 


-4.44(41) 


-4.61(61) 


-4.60(62) 


1.04(10) 



Table 3: Data on masses and couplings for the two fermion doublers at one-7r 
corner. Labels are the same as in table 1. 







l/W l7r )l 


GS(1t) 


GS(lTr) 


#i(br) 


R 2 {lit) 


a 


2.00(1) 


2.00(1) 


0.63(11) 


-0.59(11) 


0.10(3) 


0.10(3) 


b 


1.90()1 


2.09(1) 


1.19(14) 


-1.32(15) 


0.13(3) 


0.13(2) 


c 


1.65(2) 


2.30(3) 


1.13(16) 


-1.65(19) 


0.19(2) 


0.20(2) 


cc 


1.87(3) 


2.39(3) 


1.93(17) 


-2.36(21) 


0.22(2) 


0.21(2) 


d 


1.44(3) 


2.54(6) 


2.13(20) 


-3.68(35) 


0.37(4) 


0.36(4) 


D 


1.45(18) 


2.78(22) 


1.96(22) 


-3.82(57) 


0.32(4) 


0.33(4) 


e 


1.48(3) 


2.41(6) 


2.04(21) 


-3.31(30) 


0.32(4) 


0.31(3) 


E 


1.50(4) 


2.40(10) 


1.85(21) 


-3.21(53) 


0.28(3) 


0.31(5) 
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